Simultaneous photon absorption as a probe of molecular interaction and hydrogen-bond 

correlation in liquids 
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We have investigated the simultaneous absorption of near-infrared photons by pairs of neighboring molecules 
in liquid methanol. Simultaneous absorption by two OH-stretching modes is found to occur at an energy higher 
than the sum of the two absorbing modes. This frequency shift arises from interaction between the modes, and 
its value has been used to determine the average coupling between neighboring methanol molecules. We find 
a rms coupling strength of 46 ± 1 em , much larger than can be explained from transition-dipole coupling, 
suggesting that hydrogen-bond mediated interactions between neighboring molecules play an important role 
in liquid methanol. The most important aspect of simultaneous vibrational absorption is that it allows for a 
quantitative investigation of hydrogen-bond cooperativity. We derive the extent to which the hydrogen-bond 
strengths of neighboring molecules are correlated by comparing the line shape of the absorption band caused 
by simultaneous absorption with that of the fundamental transition. Surprisingly, neighboring hydrogen bonds 
in methanol are found to be strongly correlated, and from the data we obtain a hydrogen-bond correlation 
coefficient of 0.69 ±0.1 2. 



When two molecules are in close proximity, there is a fi- 
nite probability for them to simultaneously absorb a single 
photon and 'share' its energy. For vibrational transitions, this 
effect was first observed more than fifty years ago in the near- 
infrared absorption spectra of gas mixtures [1], and later in 
crystalline hydrogen chloride [2] and other condensed-phase 
systems J3. Hi- Simultaneous absorption is observed as an 
absorption peak at the sum frequency of the two vibrational 
modes, and arises from a process in which one near-infrared 
photon is absorbed, and both of the molecules involved be- 
come vibrationally excited. It is obvious that simultaneous 
photon absorption can occur only if the two molecules in- 
volved 'sense' each other's presence, that is, if they inter- 
act JH. Here, we show that simultaneous vibrational absorp- 
tion in fact provides a unique probe of intermolecular inter- 
actions in liquids. It can be used to determine the vibrational 
interactions between neighboring molecules, and more impor- 
tantly, it can be used to determine quantitatively the corre- 
lation between the strengths of neighboring hydrogen bonds. 
This hydrogen-bond correlation (or cooperativity) is believed 
to be a key property of hydrogen-bonded liquids flSSl. but 
has as yet been difficult to investigate experimentally. 

We have investigated the simultaneous absorption by neigh- 
boring CD 3 OH molecules in liquid CD3OH/CD3OD mix- 
tures. By measuring on isotopically diluted solutions of 
CD3OH in CD3OD, the concentration of OH- - OH pairs 
can be varied continuously without any change in the struc- 
tural properties of the liquid. The absorption spectra were 
measured using a Fourier-transform infrared spectrometer 
(Bruker Vertex 70) with a spectral resolution of 0.5 cm" . 
Samples were prepared by mixing appropriate amounts of 
CD3OH and CD3OD (Euriso-Top SA, >99.80%D, HDO/D 2 
content<0.03%). The weak background absorption caused 
by overtones and combination modes of CD-stretching vibra- 
tions y[] was eliminated by subtracting the spectrum of pure 
CD3OD measured in the same sample cell as used for the mix- 



tures. Samples were kept between CaF 2 windows separated 
by either a 2.2 or 10.0 mm thick spacer. 
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FIG. 1: (a) Absorption spectra of CD3OH/CD3OD mixtures with 
increasing molar fractions xcd,OH (sample-cell thickness 10.0 mm, 
CD3OD absorption subtracted), (b) Same spectra, scaled to the OH- 
stretch overtone (v = — > 2) intensity at 6300 cm~ 1 . 



Figure [T^ shows the near-infrared absorption spectrum of 
a series of CD3OH/CD3OD mixtures with increasing molar 
fraction jccd 3 oh (CD3OD background subtracted). The most 
prominent feature of these spectra is the band centered at 
~6450 cm -1 , which is caused by v = — > 2 overtone absorp- 
tion of the OH-stretch mode of CD3OH yfl. The intensity of 
this band increases linearly with the OH-stretch concentration 
■*cd 3 oh- Closer inspection shows that as the concentration of 
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OH groups increases, a second absorption feature appears at 
higher frequency. This can most clearly be seen when the 
spectra are scaled to the overtone absorption and overlayed, 
see Figure QJ). As the feature appears and grows with increas- 
ing CD3OH concentration, it can be assigned to simultaneous 
absorption by more than one CD3OH molecule QQ]. 

By subtracting the absorption spectrum at low CD3OH con- 
centration (where the simultaneous absorption is relatively 
smaller) from the spectra at higher xcd 3 oh, the simultaneous 
absorption peak can be observed separately from the more in- 
tense overtone band. The simultaneous absorption spectrum 
obtained in this way is shown in Figure|2](solid curve). Iden- 
tical simultaneous absorption spectra (apart from an overall 
scaling factor) are obtained for all values of jccd 3 oh up to 0.3. 
As the center frequency of the simultaneous absorption band 
(^6700 cm" 1 ) is close to twice the fundamental OH-stretch 
frequency of CD3OH, it seems logical to assign the peak 
to simultaneous absorption of one photon by the OH-stretch 
modes of two CD3OH molecules in close proximity |4J]. To 
verify that indeed two, and not more, CD3OH molecules are 
involved, we have determined the concentration dependence 
of the simultaneous absorption (Figure[3]). Whereas the over- 
tone absorption (solid points) varies linearly withxcD 3 OH con- 
firming that it is due to "one-molecule absorption", the simul- 
taneous absorption intensity (open points) varies quadratically 
with xcd 3 oh, which proves that the simultaneous absorption 
involves pairs of CD3OH molecules ("two-molecule absorp- 
tion"). It may be noted that the quadratic xcd 3 oh dependence 
can already be seen directly from the linear dependence of the 
simultaneous absorption feature in the scaled data of Fig.QJ. 
We conclude that the observed simultaneous absorption peak 
must be due to the absorption of photons by the OH-stretching 
modes of pairs of CD3OH molecules, since all other normal 
modes of CD3OH have frequencies too low to give a sum fre- 
quency of ~6700 cm" 1 1 9]. 




6300 6400 6500 6600 6700 6800 6900 7000 7100 7200 



Frequency V (cm ) 

FIG. 2: Solid curve: simultaneous absorption band, obtained by sub- 
tracting the absorption with xcd,OH = (scaled to overtone ab- 
sorption) from that with xcd,OH = 0.2. Dashed curve: v = — > 1 
absorption spectrum, multiplied horizontally by a factor of 2. The 
simultaneous absorption occurs at a frequency 36 cm" 1 higher than 
twice the fundamental frequency. 
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FIG. 3: Log-log plot of the intensities of the v = — > 2 overtone 
band and of the simultaneous absorption band vs. xcd 3 oh (the si- 
multaneous absorption points have been vertically shifted by 1 unit 
for better comparison). Also shown are lines with slopes of 1 and 
2, corresponding to linear and quadratic concentration dependence, 
respectively. 

As already stated, simultaneous absorption can arise only 
if there is interaction between the vibrational modes involved. 
In the present case, the interaction required is between the OH 
groups of two methanol molecules, which should therefore be 
in sufficiently close proximity. It is easily demonstrated that 
only pairs of OH groups that are direct neighbors can con- 
tribute significantly to the observed simultaneous absorption 
band. The intensity of simultaneous absorption is proportional 
to the square of the interaction between the oscillators jj§]. As 
a consequence, the simultaneous-absorption intensity due to 
dipolar interaction (the dominant interaction mechanism lllOll ) 
between two OH groups decays with their distance as r" 6 ; and 
the non-dipolar contributions to the interaction (occurring for 
instance through the OH- • • OH hydrogen bonds ED, EH) de- 
cay even faster with distance. Hence, the contribution of OH 
pairs that are next-nearest neighbors (or even farther apart) 
is negligible, and the simultaneous absorption is caused by 
neighboring OH groups only. 

Interestingly, the simultaneous absorption occurs not at ex- 
actly twice the OH-stretch frequency (2 x 3343 = 6686 cm" 1 ), 
but at a frequency which is higher (6722 cm" 1 ). This can 
be seen clearly in Figure [2] where the simultaneous band is 
compared with the fundamental absorption band multiplied 
horizontally by 2 (dashed curve). This dashed curve corre- 
sponds to the hypothetical case of simultaneous absorption 
by two non-interacting CD3OH molecules (note that because 
of the low xcd 3 oh, the fundamental and overtone bands are 
due to absorption by isolated CD3OH molecules, the con- 
centration of which is much higher than that of paired ones). 
Apparently, the interaction between the paired OH groups is 
large enough to change the energy of the simultaneously ex- 
cited state significantly. From the position of the simultane- 
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ous absorption band with respect to the doubled fundamen- 
tal frequency (which corresponds to the situation in absence 
of coupling), we can derive the (average) magnitude of the 
coupling between the OH-stretching modes. This coupling is 
dominated by the contribution lowest-order in the two OH- 
stretch displacements |9|. and we can write it as V = $qiq2, 
where q\ and q2 are the normal-coordinate displacement op- 
erators of the neighboring OH-stretch modes, and p the cou- 
pling strength. We treat the coupling as a perturbation, de- 
noting by \ij) the unperturbed state having the OH-stretch 
mode of molecule A in the v = i vibrational level, and the 
OH-stretch mode of neighboring molecule B in the v = j vi- 
brational level. Using harmonic approximations for the wave 
functions (HHQ, we have (l%i# 2 |H) = 0, so that the first- 
order correction to the energy of the simultaneously excited 
state 1 11) vanishes lfl4ll . The second-order energy correction 
with respect to the situation without coupling is given by 01511 
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where E^' is the unperturbed energy of state \ij). Because of 
the i — * i ± 1 selection rule for the operator, the V matrix 
element in the numerator is nonzero for only a few states. Of 
these, the only ones sufficiently close in energy to contribute 
significantly to the above sum are the |02) and |20) overtone 
states, see Figure[4] Hence, we have 
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where we have used harmonic vibrational wave functions to 
calculate the matrix element 11511 . and where (. . .) denotes en- 
semble averaging over all OH pairs contributing to the simul- 
taneous absorption band. Since only nearest neighbors con- 
tribute significantly, this averaging is essentially over the rel- 
ative orientation distribution of two neighboring OH groups. 

(2) 

Setting (AZsjj ) in the above expression equal to the observed 
shift of 36 ±2 cm -1 , and substituting the unperturbed values 
and E^ (see Fig. 3J, we obtain the root-mean-square 
coupling constant ((3 2 ) 1 / 2 =46± 1 cirr 1 !!^. 
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FIG. 4: Energy-level diagram (not to scale) for the OH-stretch two- 
exciton states. On the left the energy levels in absence of interaction 
between the two OH groups (isolated OH groups), on the right the 
energy levels when the molecules are neighbors. 



This experimentally determined OH-OH interaction may be 
compared to the transition-dipole coupling between the two 
OH groups, which is in general the dominant contribution to 
intermolecular vibrational couplings lllOl Il7ll . The transition- 
dipole coupling strength between two OH groups A and B 

is EH 
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where Tab is the distance vector between the transition dipoles 
fi A and n B . The values of /ja=L>b = 0.264 D and tab = 2. 84 A 
are known 111 811 . but the angle between the transition dipoles 
of the two CD3OH molecules will vary from pair to pair. To 
obtain an upper limit for the transition-dipole coupling, we 
calculate Pdipole for the case of parallel transition dipoles, for 
which we find pdipole = —30 cm" 1 . This upper limit has 
a magnitude much smaller than the observed rms value of 
46 cm" 1 . Hence, the coupling between two neighboring OH- 
groups is much larger than can be explained from transition- 
dipole coupling alone, and apparently other coupling mech- 
anisms, like hydrogen-bond mediated coupling 111 li Il2l 1 1711 . 



contribute strongly to the molecular interaction in methanol. 

Perhaps the most important and surprising aspect of si- 
multaneous vibrational absorption is that it can be used to 
study correlation among hydrogen bonds. There is strong 
evidence that in hydrogen-bonded liquids, hydro- 

gen bonding occurs in a cooperative way: an OH group 
acting as an acceptor in a strong hydrogen bond in turn 
tends to form strong hydrogen bonds as a donor, and vice 
versa. It is well known that the OH-stretch frequency of a 
hydrogen-bonded OH group is proportional with the strength 
of the hydrogen bond (low OH-stretch frequency correspond- 
ing to short OH- • O hydrogen-bond distance) lfl9l, I20I1 . In 
hydrogen-bonded liquids, there generally exists a distribu- 
tion of hydrogen-bond strengths, leading to an OH-stretch 
absorption band that is inhomogeneously broadened 1I20I 
I2III . its width reflecting the width of the distribution of 
hydrogen-bond lengths. Because of the proportionality be- 
tween hydrogen-bond length and OH-stretch frequency, a cor- 
relation between neighboring hydrogen bonds should give rise 
to an identical correlation between the OH-stretch frequencies 
of neighboring OH groups. Hence, by determining the corre- 
lation between the OH-stretch frequencies Va and Vs of two 
neighboring OH groups, we can directly determine the cor- 
relation between their hydrogen bonds. To what extent two 
frequencies Va and Vb are correlated can be quantitatively de- 
termined by comparing the width of the distribution of the 
sum Va + Yb with that of the fundamental distributions of Va 
and Vb [22], that is, by comparing the widths of the simulta- 
neous and fundamental absorption bands. We define the cor- 
relation coefficient between the frequencies Va and Vb in the 
usual manner as 



((va-(v»(vb-(v))> 



(4) 



where (. . .) denotes ensemble averaging, and where we have 
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defined (v) = (Va) = (vb) and o v = Va = Vs (note that this 
is possible because the probability distributions of Va and Vb 
are identical). It is easily shown that the correlation coefficient 
is related to the ratio of the widths of the fundamental and 
simultaneous absorption bands by [22] 

^™ = V2 + 2^, (5) 

a result that is independent of the functional form of the Va,b 
probability distribution. The familiar limiting cases of this 
equation are: completely independent frequencies (r = 0), 
where ol A+VB = g Va + o Vg = 2o^ so that o Va +v b = V2a v , 
and completely correlated frequencies (r = 1), where Va = Vb 
so that Gv A +v B — 02v A — 2c v . The simultaneous absorption 
bands predicted for these two limiting r values are drawn as 
dotted curves in Figure|5] The observed simultaneous absorp- 
tion band lies between these limiting cases, which implies 
that the OH-stretching frequencies of neighboring methanol 
molecules are partially, but not completely correlated. From 
the observed ratio o Va+Vb /o v of the widths of the simultane- 
ous and fundamental absorption bands we can determine r by 
means of Eq. ©. We observe o Va + Vb /o v — 1.84±0.06, from 
which a correlation coefficient r = 0.69 ± 0. 12 is obtained. To 
our knowledge, this constitutes the first experimental determi- 
nation of the correlation coefficient between neighboring hy- 
drogen bonds in a liquid, a fundamental parameter that char- 
acterizes hydrogen-bond cooperativity in a quantitative way. 
As yet, no theoretical prediction is available for this number, 
and we hope the results presented here will stimulate work in 
this direction. 
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FIG. 5: Solid curve: simultaneous absorption band. Dashed curve: 
fundamental absorption band. Both bands have been centered at the 
origin. Dotted curves: fundamental absorption band, horizontally 
multiplied by a factor of \fl and by a factor 2. These two curves 
are the simultaneous absorption bands expected if the two coupled 
OH-oscillators are completely uncorrelated (r = 0), and completely 
correlated (r = 1), respectively. The experimental simultaneous ab- 
sorption band is between the r = and r = 1 curves, indicating that 
the OH-frequencies are partially correlated. 
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